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This  paper  presents  a  grid-connected  photovoltaic  system  (PV)  used  as  a 
shunt  active  power  filter  (SAPF)  to  provide  the  power  factor  correction, 
harmonic  elimination,  reactive  power  compensation  and  to  simultaneously 
supply  power  from  a  PV  system  to  the  utility.  A  direct  power  control  (DPC) 
method  is  used  for  controlling  the  system  to  feed  the  photovoltaic  energy  in 
synchronization  with  grid  and  provide  power  quality  improvement.  The  PI 
parameters  of  DC-link  voltage  controller  are  tuned  using  the  Particle  Swarm 
Optimization  (PSO)  algorithm  without  the  need  for  an  exact  mathematical 
model  of  system.  This  PI-PSO  controller  gives  better  results  for  robustness, 
harmonic  minimization  and  reduces  the  overshoot  and  undershoots  of  PI 
controller.  The  overall  control  of  system  is  tested  in  Matlab/Simulink 
environment.  Then,  the  simulations  results  demonstrate  the  robustness  and 
feasibility  of  proposed  method. 
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1.  INTRODUCTION 

The  solar  energy  is  the  most  popular  source  among  the  alternative  energy  resource  without  emission 
of  pollutants  energy  conversion  is  done  [1],  In  an  effort  to  use  the  solar  energy  effectively,  a  great  deal  of 
research  has  been  done  on  the  grid-connected  photovoltaic  generation  systems.  In  PV  systems  connected  to 
the  grid,  the  inverter  that  converts  the  output  direct  current  (DC)  of  the  solar  modules  to  the  alternate  current 
(AC)  is  receiving  increased  interest  in  order  to  generate  power  to  utility  [2], 

Today,  the  widespread  use  of  nonlinear  loads  in  power  systems  causes  circulation  of  harmonics  and 
reactive  currents  in  the  distribution  lines  [3]  [4].  If  the  grid -connected  PV  system  is  applied  to  nonlinear 
loads,  the  power  quality  is  relatively  poorer.  To  avoid  these  problems,  the  grid-connected  PV  system  should 
not  only  supply  active  power  to  the  system  via  MPPT,  but  also  improve  the  power  quality  (low  THD  and 
unity  power  factor)  [5], 

By  modifying  the  control  of  the  inverter,  active  filtering  feature  can  also  be  included  in  the  PV 
interfacing  inverter  and  thus  saves  the  cost  of  additional  active  filter.  Various  control  types  of  shunt  active 
filter  were  presented  in  the  literature  [6].  In  this  article,  we  propose  the  direct  power  control  (DPC).  The  DPC 
strategy  is  based  on  the  instantaneous  active  and  reactive  power  comparators.  The  active  power  command  is 
provided  from  a  dc-bus  voltage  controller  block,  while  the  reactive  power  command  is  directly  given  from 
the  outside  of  the  controller.  Errors  between  the  commands  and  the  estimated  feedback  power  are  input  to  the 
hysteresis  comparators  [7]  [8]. 

In  the  DPC  technique,  the  active  power  is  delivered  from  the  outer  proportional-  integral  (PI) 


Journal  homepage:  http://iaescore.com/journals/index.php/IJPEDS 


Int  J  Pow  Elec  &  Dri  Syst 


ISSN:  2088-8694 


□  611 


controller.  This  PI  controller  is  formulated  as  an  optimization  problem[9].  In  this  work,  we  propose  the  use  of 
the  particle  swarm  optimization  (PSO)  method  to  search  for  the  optimal  parameters  of  a  PI  controller.  DPC 
based  on  classical  PI  regulator  and  PSO  algorithm  is  developed  and  a  comparison  study  is  presented.  To 
show  the  performances,  a  numerical  simulation  is  performed. 


2.  SYSTEMS  CONFIGURATIONS 

The  structure  of  the  grid-connected  PV  system  is  shown  in  Figure.  1.  It  consists  of  a  photovoltaic 
array  connected  through  a  dc-dc  boost  converter  to  a  three-phase  inverter  that  is  connected  to  a  grid  through  a 
simple  filter  and  nonlinear  load. 


Figure  1 .  Grid-Connected  Inverter  Interfaced  with  PV  System 


2.1.  PV  Array  Model 

The  PV  module  is  the  unity  of  base  for  the  construction  of  a  P  V  array.  We  have  used  the  PV  module 
model  of  a  single  diode.  This  model  is  the  most  commonly  used  due  to  good  results  [10].  The  current-voltage 
relationship  is  given  by: 


PV 
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Where: 

Iph:  the  current  generated  by  the  incident  light;  Is:  the  reverse  saturation  or  leakage  current  of  the  diode;  Rs: 
the  equivalent  series  resistance;  Rsh:  the  equivalent  parallel  resistance;  q:  the  electron  charge  (1.60217646  x 
10-19  C);  k:  the  Boltzmann  constant  (1.3806503  x  10-23  J/K);  A:  Quality  Factor  of  diode. 


Ipy 


Figure  2.  Equivalent  Circuit  of  a  PV  Cell 
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Figure  3.  (a)  I-V  Characteristics,  (b)  P-V  Characteristics 


The  PV  model  is  simulated  using  the  Sun  Power  SPR-305E-WHT-D  PV  panel.  The  simulated  I-V 
and  P-V  characteristics  of  the  SPR-305E-WHT-D  PV  module  at  standard  climatic  conditions  are  shown  in 
Figure.  3.  (a)  and  (b),  respectively. 

2.2.  Boost  Converter  Modeling 

Figure.4  shows  the  equivalent  circuit  of  the  DC-DC  boost  converter.  These  converters  can  be 
described  using  ordinary  differentiation  equations  as  follows: 


L 


Figure  4.  Control  Scheme  of  DC-DC  Boost  Converter 
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D:  is  the  duty  cycle. 

The  switch  is  on  for  a  period  equal  to  DT,  and  the  switch  is  off  for  an  interval  time  equal  to  ( 1  -  D)  T. 

2.3.  Maximum  Power  Point  Tracking 

Several  publications  on  MPP  methods  are  regularly  published  in  the  literature.  For  our  study,  we 
chose  the  MPPT  algorithm  of  P&O.  P&O  method  is  one  of  the  popular  methods  to  track  the  maximum-power 
point  and  probably  the  most  frequently  used  in  practice,  mainly  due  to  its  easy  implementation.  [11]  [12] .  The 
working  principle  of  P&O  is  depicted  by  the  flow  chart  in  Figure.  5. 
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Figure  5.  Flowchart  of  the  P&O  Algorithm 


2.4.  Active  Power  Filter 

The  structure  of  the  active  filter  consists  essentially  of  two  parts,  a  power  and  a  control  part.  The 
power  part  consists  of  an  inverter,  a  coupling  filter  and  a  source  of  energy.  The  control  part  is  used  to  control 
the  switching  of  semiconductor  elements  forming  the  inverter  from  the  power  part.  Au  through  appropriate 
control  strategies,  it  is  possible  to  generate  harmonic  signals  at  the  output  of  the  inverter  to  offset  those 
present  on  the  power  grid,  such  as.  [13] 

4  =*/+*/  (6) 

is:  current  of  source.;  z, :  current  of  the  load;  if.  current  compensation. 


3.  DIRECT  POWER  CONTROL 

The  main  idea  of  the  direct  power  control  (DPC)  originally  proposed  by  Ohnishi  (1991)  and  further 
developed  by  Noguchi  and  Takahachi  in  1998,  is  similar  to  the  direct  torque  control  (DTC)  of  induction 
machines.  Instead  of  flux  and  torque,  the  active  power  (P)  and  reactive  power  (q)  are  selected  as  two 
instantaneous  magnitudes  to  control  [14]  [7].  The  figure  .2  shows  the  configuration  of  the  direct  power 
control  without  voltage  sensor  for  a  three  phase  inverter. 


Figure  6.  General  Structure  of  Shunt  Active  Filter  Controlled  by  DPC 


The  estimated  values  of  (ps)  and  (qs)  in  terms  of  the  three  phase  line  currents  and  voltages  can  be 
derived  as:  [14] 

Ss=Ps+jQs  (7) 

^5  ~^sa'^sa  "^^sb’hb  ^ sc sc 
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Qs  vsb  A,  )  xsa  +  ( vsc  v sa )  .isj}  +  ( vsa  vsb )  ,isc  J 


(9) 


The  phase  of  the  power-source  voltage  vector  is  converted  to  the  digitized  signal  0N  Error! 
Reference  source  not  found..  For  this  purpose,  the  working  plane  (a,  P)  is  divided  into  12  sectors  in  Figure. 
7.  The  sectors  can  be  numerically  determined  by  the  following  equation.  [7]. 

Or:  N=  1,2, 3....  12  is  the  number  of  the  sector 


6  =  Arctg 
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Figure  7.  Plant  (a,  P)  Divided  in  12  Sector 


3.1.  Hysteresis  control  and  switching  table 

The  references  of  reactive  power  (qref)  and  active  power  (pref)  are  compared  with  the  estimated 
reactive  (qs)  and  active  ( ps )  values.  Errors  between  the  commands  and  the  estimated  feedback  power  are 
input  to  the  hysteresis  comparators  and  digitized  to  the  signals  Sp  and  Sq.  The  output  signal  of  the  active 
power  controller  is  defined: 

(Pref-Ps  <~hp=>Sp=0  (n) 

|  Pref  ~  Ps  >  ~hp  =>  Sp  =  1 

And  similarly  for  the  reactive  power  controller  as: 

\qref-qs<-hq  =>Sq  =0 

[  (Inf  -‘Is  >~hq  =>Sq  =1 

For  all  sectors  the  switching  table  proposed  is  represented  in  the  Table  1. 


Table  1.  All  Sectors  Switching  Table 
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3.2.  PI  controller 

The  regulation  of  the  continuous  voltage  of  inverter  is  ensured  by  a  regulator  of  PI  type.  The  error 
between  the  voltage  of  the  capacitor  Vdc  and  the  reference  voltage  Vdcief  is  used  as  an  input  to  PI  Controller 
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and  the  output  is  the  reference  active  power  pref.  The  transfer  function  of  the  PI  Controller  is: 

G(s)  =  kp  +  —  (13) 

Where,  kP  is  the  proportional  constant  gain  and  k,  is  the  integral  constant  gain.  These  gains  are  tuned 
by  using  Particle  Swarm  Optimization  (PSO). 


4.  PARTICLE  SWARM  OPTIMIZATION 

The  particle  swarm  optimization  (PSO)  algorithm  was  first  described  in  1995  by  Kennedy  and 
Eberhart  [15]  [16].  The  technique  has  evolved  greatly  since  then.  PSO  is  a  stochastic,  population  based 
evolutionary  computer  algorithm  for  problem  solving.  In  a  PSO  system,  a  swarm  of  individuals  (called 
particles)  fly  through  the  search  space.  Each  particle  represents  a  candidate  solution  to  the  optimization 
problem.  The  position  of  a  particle  is  influenced  by  the  best  position  visited  by  itself  (i.e.,  its  own  experience) 
and  the  position  of  the  best  particle  in  its  neighbourhood  (i.e.,  the  experience  of  neighboring  particles).  When 
the  neighbourhood  of  a  particle  is  the  entire  swarm,  the  best  position  in  the  neighbourhood  is  referred  to  as 
the  global  best  particle,  and  the  resulting  algorithm  is  referred  to  a  global  best  PSO.  When  smaller 
neighbourhoods  are  used,  the  algorithm  is  generally  referred  to  a  local  best  PSO. 

The  performance  of  each  particle  (i.e.,  how  close  the  particle  is  from  the  global  optimum)  is 
measured  using  a  fitness  function  that  varies  depending  on  the  optimization  problem.  Each  particle  in  the 
swarm  is  represented  by  the  following  characteristics  [15]  [17]: 

x,  is  the  current  position  of  the  particle;  v,  is  the  current  velocity  of  the  particle;  y,  is  the  personal  best  position 
of  the  particle  is  the  neighbourhood  best  position  of  the  particle.  The  personal  best  position  of  particle  i  is  the 
best  position  visited  by  particle  i  so  far.  Let  F  denote  the  objective  function.  Then  the  personal  best  of  a 
particle  at  time  step  t  is  updated  as: 

f  y  ■  (t)  if  F{xAt  +  \)>F{yAt)) 

v.(f  +  l)  =  J  ■r,v  (14) 

[Xj (t  + 1) if  F(Xj(t  +  Y)  <  F (yi (f )) 

For  the  gbest  model,  the  best  particle  is  determined  from  the  entire  swarm  by  selecting  the  best 
personal  position.  If  the  position  of  the  global  best  particle  is  denoted  by  the  vector  then: 

Where: 


}  =  min{F(y0(t)), . ,F  (ys(t))}  (15) 

And:  s  denotes  the  size  of  the  swarm. 

The  velocity  update  step  is  specified  for  each  dimension 
Z  :  j  ^,{1, . ,Nrf} 

Hence,  v,  ■  represents  the  jth  element  of  the  velocity  vector  of  the  ith  particle.  Thus  the  velocity  of 
particle  i  is  updated  using  the  following  equation: 

\’i  j  =  co.Vj  j  (f)  +  c,  .A,  +  c2 ,A2  ( 1 6) 

Where: 

Ar  =^.(^.(0-^(0)  (17) 

a2  =  r2j-(y/(0— *ij(0)  (is) 

co  is  the  inertia  weight,  q  and  c2  are  the  acceleration  constants,  and  random  coefficients  distributed  as: 

\j  and  r2  ]  e{0.l} 

The  position  of  particle  i,  x,  is  then  updated  using  the  following  equation: 
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xi(t  +  ])  =  xl(t)  +  vl(t  +  ])  (19) 

This  process  is  repeated  until  a  specified  number  of  iterations  is  exceeded,  or  velocity  updates  are 
close  to  zero.  The  quality  of  particles  is  measured  using  a  fitness  function  which  reflects  the  optimality  of  a 
particular  solution.  The  following  steps  summarize  the  basic  PSO  algorithm: 

4.1.  PSO  Algorithm 

For  each  particle  i=l, . s  do 

Randomly  initialize  x, 

Randomly  initialize  v,  (or  just  set  v,  to  zero ) 

Set  yi=Xj 
endfor 
Repeat 

For  each  particle  i=l,  ...,s  do 
Evaluate  the  fitness  of  particle  i,  f(xj 
Update  y ,  using  equation  (16) 

Update  {’  using  equation  (18) 

For  each  dimension  j=  l,...,Nd  do 
Apply  velocity  update  using  equation  (19) 

Endloop 

Apply  position  update  using  equation  (21 ) 

Endloop 

Unit l  some  convergence  criteria  is  satisfied 

4.2.  Objective  function  formulation 

In  this  study,  the  PSO  algorithm  is  proposed  to  tune  the  PI  parameters  of  voltage  controller  by 
minimizing  the  cost  function.  So  that,  the  proportional  (kp)  and  integral  (ki)  gains  of  PI  controller  is  defined 
as  the  problem  parameters  that  should  be  optimized. 

Thus,  each  individual  of  population  is  composed  of  the  controller  parameters  kp  and  ki.  Then,  the 
optimization  problem  of  PI  controller  is  formulated  as  the  minimization  of  the  integral  time  absolute  error 
ITAE  and  of  the  maximum  overshoot.  Therefore,  the  objective  function  is  defined  as: 


l  =  U;^rrt[dt 

t= l 

k  <  k  <  k  k  <  k  <  k 

'v/?mm  —  —  A'z? max  ’  'vimin  —  /  —  m 


Where: 

n:  is  the  total  number  of  points  for  which  the  optimization  is  carried  out. 
ts :  is  the  time  rang  of  simulation  or  simulation  period. 


(20) 


Figure  8.  Insertion  of  PSO  Approach  with  the  PI  Controller 


5.  SIMULATION  RESULTS  AND  ANALYSIS 

The  system  depicted  in  Figure.  1  and  6  has  been  simulated  using  the  SimPowerSystem  with  S- 
function  of  Matlab/Simulink  installed  on  a  computer  system  with  Intel  Core  i5-4210  CPU,  1.70  GHz  speed 
and  4GO  of  memory.  Parameters  used  in  simulation  are  as  follows: 
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Table  2.  Parameters  of  PSO  Algorithm 


Parameter 

Values 

Swarm  size 

10 

Number  of  iterations 

100 

rl 

0.3 

r2 

0.95 

cl 

1.45 

c2 

1.45 

Table  3.  Parameters  of  DC-DC  Boost  Converter 

Parameter 

Values 

/switching 

5kHz 

Output  voltage 

400V 

Inductor  L 

10  mH 

Capacitor  C 

lOOpF 

Table  4.  SAPF  Parameters. 

Parameter 

Values 

Source  voltage 

125  V 

Supply  frequency 

50Hz 

Source  impedance  Rs,  Ls 

O.in,  0.1  mH 

Filter  impedance  Rf,  Lt 

0.1  a,  2.6mH 

Line  impedance  Ri,  Li 

0.1  a,  0.5mH 

DC-link  capatitor  C 

2.2pH 

Diode  rectifier  load  Rl,  Ll 

20Q,  2mH 

The  source  current  and  its  spectrum  before  connecting  SAPF  feed  by  PV  system  are  illustrated  in 
Figure  9.  The  results  show  that  the  system  is  affected  by  the  nonlinear  load  during  the  operation.  Where,  the 
Total  Harmonic  Distortion  (THD)  of  the  supply  current  is  high  and  equals  27.57%. 
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J _ 1  ■ 
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(b) 


Figure  9.  (a)  Supply  Current  before  Harmonics  Compensation,  (b)  Harmonic  Spectrum. 


To  verify  the  performances  and  robustness  of  the  proposed  control  and  the  research  algorithm  of 
maximum  power  point  tracking  (MPPT),  a  variable  profile  irradiance  to  sweep  all  modes  of  operation  of  our 
system  is  done.  Figure  10  shows  the  solar  irradiance  level  starts  from  0  W/m2,  then  increases  to  600  W/m2, 
after  that  decreases  to  200  W/m2  and  finally  increases  to  1000  W/m2. 
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(a) 


(a) 


(b) 


Figure  10.  (a)  Solar  Radiation,  (b)  Current  PV  Array. 


Figure  11.  Simulation  Results  of  the  Proposed  DPC  with  PI  Controller:  (a)  Source  Current,  (b)  Source 
Voltage  and  Source  Current,  (c)  Filter  Current,  (d)  Reactive  Power,  (e)  DC-Link  Voltage  vdc,  (f)  Active 

Power  and  PV  Power. 
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Figure  12.  Simulation  Results  of  the  Proposed  DPC  with  PI-PSO  Controller:  (a)  Source  Current,  (b)  Source 
Voltage  and  Source  Current,  (c)  Filter  Current,  (d)  Reactive  Power,  (e)  DC-Link  Voltage  Vdc,  (f)  Active 

Power  and  PV  Power. 


Figure  11  and  12  shows  the  system  with  PI  and  PI-PSO  controllers.  The  source  current  (Is),  source 
voltage  (Vs),  filter  current  (If),  DC  bus  voltage  (Vdc),  reactive  power  (qs),  active  power  (Ps)  and  PV  power 
(Ppv)  are  depicted.  The  results  show  that  the  the  network  current  become  perfectly  sinusoidal  after  the 
control  application  and  remain  in  phase  with  the  corresponding  voltages. The  change  in  radiation  don’t  affect 
DC  link  voltage  but  change  the  line  side  active  power  and  the  amplitude  of  the  line  current.  The  reactive 
power  is  always  zero. 

We  can  observed  from  [0.3s-0.5s]  and  [0.8s- Is]  ,  the  opposite  phase  between  voltage  and  current 
source  and  the  negative  sign  of  the  utility  active  power  (Ps),  meaning  the  filter  current  has  information  of 
harmonic  and  PV  array  currents  to  ensure  elimination  of  harmonic  and  injection  of  current  to  the  load.  The  dc 
link  voltage  returns  to  its  reference  value  in  few  milliseconds,  without  overshoot  and  undershoots,  as  shown 
in  Figure  12e,  demonstrates  the  superiority  of  the  PI-PSO  controller  compared  to  its  counterpart  PI  controller. 

Figure  13  show  harmonic  spectrum  of  source  current  at  0  W/m2.  The  active  filter  decreases  the  total 
harmonic  distortion  (THD)  from  27.57%  to  1.83%  with  PI  controller.  However,  with  PI-PSO  controller  the 
THD  is  increased  to  1.28%  which  proves  the  effectiveness  of  the  proposed  controller.  Performance  of 
proposed  control  is  evaluated  in  term  of  THDi  level  behavior  under  fast  changing  irradiation  as  well  as 
changing  operation  modes  as  shown  in  Table  5. 
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(a)  Fundamental  (50Hz)  =  11.1  ,  THD=  1.28%,  (0  W/m2) 


(b)  Fundamental  (50Hz)  =  11.1  ,  THD=  1.83%,  (0  W/m2) 
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Figure  13.  Harmonic  Spectrum,  (a)  With  PI-PSO  Controller,  (b)  With  PI  Controller. 


Figure  14.  Error  Function  Evolution  with  Iterations  for  PSO  Algorithm. 


Table  5.  Comparisons  of  DPC  Indexes  between  Used  and  Unused  PSO  Algorithm 


kp 

ki 

THD  (0  W/m2) 

THD  (600  W/m2) 

THD  (200  W/m2) 

THD  (1000  W/m2) 

PI  0.2 

75.63 

1.83% 

6.17% 

3.21% 

1.66% 

PI-PSO  0.365 

23.547 

1.28% 

4.25% 

2.18% 

1.16% 

The  simulation  result  shows  that  the  improved  particle  swarm  optimization  PI  control  has  a  smaller 
overshoot  during  radiation  change,  low  rise  time  and  low  source  current  THD  compared  to  its  counterpart  PI 
controller.  It  can  be  concluded  that  the  DPC  with  PI-PSO  gives  betters  results  for  robustness,  harmonic 
minimization  and  inject  the  maximum  active  power  available  from  the  PV  array  into  the  load  and/or  grid. 


6.  CONCLUSION 

In  this  paper,  a  PSO  based  direct  power  control  for  a  multifunctional  grid  connected  photovoltaic 
system  is  investigated,  the  simulation  results  show  a  significant  enhancement  by  applying  the  DPC  with  PI- 
PSO  controller  in  comparison  with  classical  DPC  in  termes  of  response  time,  overshoot,  THD  and  stability 
the  systeme  under  irradiation  changes. The  optimization  of  the  gains  of  PI  controller  for  dc  link  voltage  by 
using  PSO  algorithm  presents  a  perfect  efficiency  to  reduce  the  overshoot  and  undershoot  of  PI  controller  and 
gives  better  results  for  robustness  and  harmonic  minimization. 
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